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INTRODUCTION 

\ 

•  In  a  free  electron  laser  a  relativistic  electron  beam  is 
bunched  under  the  action  of  the  ponderomotive  potential  and 
is  forced  to  radiate  In  close  phase  with  the  input  wave. 
Until  recently,  most  theories  of  the  FEL  have  dealt  solely 
with  electron  beams  of  Infinite  transverse  dimension 
radiating  only  one-dimensional  E.M.  waves  (plane  waves).1 
Although  these  theories  describe  accurately  the  dynamics  of 
the  electrons  during  the  FEL  interaction  process,  neither  the 
three  dimensional  nature  of  the  radiated  fields  nor  Its 
non-monochomatlc  features  can  be  properly  studied  by  them. 
As  a  result  of  this,  very  Important  practical  Issues  such  as 
the  gain  per  gausslan-spherlcal  optical  mode  In  a  free 
electron  laser  have  not  been  well  addressed,  except  through 
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a  one  dimensional  field  model  in  which  a  filling  factor 
describes  crudely  the  coupling  of  the  FEL  Induced  field  to 
the  Input  field.  *L _ ^ 

In  an  effort  to  alleviate  this  problem  a  computer  code 

Is  presently  under  development  at  UCSB  to  study  the 

three-dimentlonal  characteristics  of  fields  radiated  by  FELs. 

It  Is  based  on  the  coherent  superposition  of  electromagnetic 

pulses  radiated  by  Individual  electrons  In  the  beam.  The 

three  dimensional  fields  of  each  pulse  of  radiation  Is 

calculated  numerically  from  the  Llenard-Wlechert  solutions  to 

the  wave  equation.  Included  in  these  calculations  are  the 

effects  of  the  FEL  ponderomotlve  force  on  the  motion  of  the 

electron.  Other  approaches  to  solve  the  3D  FEL  problem  are 

2  3  4  5 

also  presently  under  development.’  *  ’ 

In  this  paper,  some  preliminary  numerical  calculations 
of  the  angular  and  spectral  characteristics  of  the  radiated 
field  of  a  constant  period  FEL  Is  presented.  Although  the 
work  is  mostly  numerical  In  nature,  the  approach  used  (l.e. 
single  particle  radiation  fields)  gives  good  understanding  of 
the  physics  of  the  FEL  wave-beam  Interaction  process. 

The  most  notable  results  presented  here  can  be 
sunmarlzed  as  follows:  Although  Individual  electron 

radiation  pulses  carry  Information  of  relatively  large 
spectral  bandwidth  and  relatively  large  angular  distribution 


patterns,  because  of  the  bunching  of  the  electron  bean  by  the 
ponderomotlve  forces,  the  coherent  superposition  of  radiation 
pulses  from  all  electrons  In  the  beam  produce,  as  a  result  of 
wave  Interference  effects,  a  narrower  cone  of  radiation  and  a 
more  monochromatic  radiation  pulse  than  the  corresponding 
effects  expected  from  a  single  electron.  It  Is  found  that 
FELs  can  produce  a  pulse  of  radiation  whose  half-angle  at 
half- power  Is  approximately  given  by  1/N,  where  N  Is  the 
number  of  periods  In  the  magnetic  wlggler  and  whose  bandwidth 
(for  long  electron  pulse  lengths)  is  approximately  given  by 
1/At,  where  At  Is  the  duration  of  the  electron  pulses. 


THE  LIENARD-WIECHERT  FIELDS  AND  EQUATION  OF  MOTION 

The  retarded  solution  of  the  wave  equation  for  an 
accelerated  point-charge  particle  moving  with  relativistic 
energy  7 me2  Is  given  by  the  well  know  Lienard-Wlechert  fields6 
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where  the  fields  are  evaluated  at  the  observation  position  x 
at  the  time  t.  The  quantities  R,  0  and  Tare  respectively 
the  radius  vector  from  the  particle  to  the  observation  point, 
the  Instantaneous  normalized  velocity  vector  and  acceleration 
vector  of  the  particle.  All  quantities  Inside  the  brackets 
must  be  evaluated  at  the  retarded  time  t'  defined  by  the 
following  relation 

t  .  t.  (3) 
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The  first  term  In  equation  (1)  describes  the  Coulomb  field  of 
a  moving  particle.  The  second  term  describes  its  radiation 
field.  The  motion  of  the  particle  Is  governed  by  the  Lorentz 
force  equation 
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It  can  be  shown8  that  the  transverse  motions  of  the  electron 
in  a  free-el ectron  laser  Is  determined  with  good  accuracy  by 
the  wlggler  field  alone.  The  longitudinal  motion  of  the 
electron,  on  the  other  hand,  is  governed  by  the  ponderomotlve 
force.  For  a  helical  wlggler  and  a  circularly  polarized 


Input  plane  wave  the  force  is  given  by: 
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where  K  *  qAoBw/2ir  me ,  and  AQis  the  period  of  the  wiggler. 
E$is  the  amplitude  of  the  input  plane  wave.  Equations(5)  and 
(6)  can  be  solved  numerically  for  the  position  z  of  the 
particle.  Using  this  result  and  equation  (3),  E  and  B  can  be 
then  evaluated  at  IT  at  time  t. 


RADIATION  FROM  A  SINGLE  ELECTRON 

The  three  dimensional  radiation  fields  from  a  single 
electron  moving  solely  under  the  influence  of  a  constant 
period  wiggler  has  been  studied  in  detail. ^  Radiation 
reaction  forces  are  small  and  do  not  affect  the  motion  of  the 
electron  as  it  radiates.  At  an  angle  of  observation  e  with 
respect  to  the  axis  of  motion  an  observer  far  away  from  the 
wiggler  will  see  coming  a  pulse  of  radiation  N  periods  long, 
of  approximate  constant  amplitude  and  phase.  The  wavelength 


of  the  radiation  field  for  a  helical  wlggler  is  given  by 
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For  small  wiggler  fields  the  radiation  can  be  shown  to  have 
an  angular  half-angle  width  9S*1A  approximately.  Also,  it 
can  be  shown  that  again  for  small  wiggler  field  most  of  the 
energy  is  radiated  within  a  narrow  band  of  wavelength  given 
very  nearly  by 
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where  N  is  the  number  of  periods  in  the  wiggler.  A  typical 
radiation  pulse  from  an  electron  moving  In  a  magnetic  wiggler 
Is  shown  by  the  broken  line  in  figure  1.  Note  that  amplitude 
and  phase  remains  very  nearly  constant  with  time.  The 
amplitude  of  the  radiation  field  is  proportional  to  the 
maximum  transverse  accelerations  given  to  the  electron  by  the 
magnetic  wiggler.  Since  the  maximum  acceleration  of  the 
electron  does  not  vary  with  positions  along  the  magnet,  then 
the  radiated  pulse  amplitude  will  not  change  with  time. 
Also,  the  phase  of  the  wave  remains  constant  because  the 
wiggler  periodicity  Is  constant  and  the  average  longitudinal 
velocity  of  the  electron  remains  constant  (i.e.  there  Is  no 


ponderomotlve  force). 
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Figure  1.  The  broken  line  represents  a  typical  radiation  field 
emitted  by  an  electron  in  a  wiggler.  The  solid  line 
represents  the  Input  plane  wave.  Both  fields  were 
observed  at  a  fixed  point  in  space. 


In  an  FEL  the  radiation  pulse  emitted  by  an  electron  Is 
different  from  that  of  a  single  electron  In  a  wiggler  alone. 
As  will  be  shown  below  the  ponderomotlve  force  affects  mainly 
the  phase  of  the  radiated  field.  The  geometrical  arrangement 
Is  illustrated  In  figure  2. 


Figure  2.  Geometrical  arrangement  corresponding  to  an  FEL 
amplifier  configuration. 


It  can  be  easily  shown  that  for  practical  FELs  the 
longitudinal  acceleration  of  the  electron  under  the  influence 
of  the  ponderomotive  force  is  much  smaller  than  the 

transverse  acceleration  of  the  electron  provided  by  the 
magnetic  wiggler. 
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The  above  ratio  can  be  interpreted  as  the  ratio  of  the  energy 
radiated  by  the  electron  in  one  optical  wavelength  to  the 
total  energy  of  the  electron.  Even  at  an  optical  power 
density  where  gain- saturation  begin  to  occur  the  above 
relation  gives 
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For  example  if  N=20  and  y-7  then  the  above  ratio  yields 
1 8z  1/ 1 8j  ssiix  to'6.  The  above  values  of  rand  N  are  the 
ones  used  to  obtain  the  results  presented  in  this  paper.  As 
a  result  of  the  above  calculation,  the  longitudinal 
acceleration  term  produced  by  the  ponderomotive  force  in  a 
FEL  does  not  contribute  significantly  to  the  amplitude  of  the 
radiated  field  and  has  thus  not  been  included  in  these 
calculations.  However,  as  it  is  indicated  in  figure  1  by  the 


broken  line,  the  effect  of  the  ponderomotlve  force  Is  to 
change  the  longitudinal  position  of  the  electron  and  thus 
change  the  phase  of  the  radiated  wave  as  shown  In  the  figure. 
In  other  words,  the  ponderomotlve  force  affects  equation  (3) 
to  a  much  greater  extent  than  It  affects  equation  (1). 


RADIATION  FROM  A  BEAM  OF  ELECTRONS  IN  A  HIGGLER 

If  there  Is  no  ponderomotlve  force,  the  radiation  fields 
produced  by  a  beam  of  electrons  moving  along  a  magnetic 
wlggler  will  depend  on  the  way  the  charge  is  distributed 
along  the  beam.  For  example,  if  the  electron  beam  contains  a 
uniform  longitudinal  charge  distribution  then,  due  to  the 
uniform  distribution  of  phases  of  the  individually  radiated 
pulses,  the  resultant  radiation  fields  will  be  exactly  zero. 
For  this  situation  only  static  electromagnetic  fields  can  be 
present.  If,  on  the  other  hand,  the  longitudinal  charge 
distribution  in  the  beam  is  random,  then  the  distribution  in 
phases  of  the  individual  radiation  pulses  will  be  also 
random.  As  a  result  of  this,  there  will  be  a  finite 
resultant  fluctuating  radiation  field  with  mean  amplitude 
proportional  where  Ne  is  the  number  of  electrons  that 

radiate  pulses  which  can  interfere  simultaneously  at  a  fixed 
point  in  space.  The  radiation  fields  produced  will  have  on 


the  average  the  angular  and  spatial  radiation  characteristics 
of  that  of  a  single  electron  as  discussed  in  the  previous 
section. 

If  the  ponderomotive  FEL  force  is  present,  then  the 
radiation  characteristics  of  the  electron  beam  can  be 

drastically  different  from  that  of  a  single  particle.  The 
amplitude  of  the  resultant  radiated  field  can  be  forced  to  be 
proportional  to  Ne,  rather  than  to  yJHe  for  random 
distribution  of  phases,  as  a  result  of  the  bunching  effects 
produced  by  the  ponderomotive  force.  This  is  the  basic 
mechanism  by  which  free  electron  lasers  work. 

RESULTS  OBTAINED 

In  this  section,  a  summary  is  presented  of  the  numerical 
analysis  of  the  three  dimensional  fields  radiated  by  a 
constant  period  free-electron  laser.  In  figure  2  is 
illustrated  the  geometrical  arrangement.  L  is  the  electron 
pulse  length  of  the  filamentary  beam  shown  as  a  heavy  dark 
line.  \Q  Is  the  magnetic  helical  wiggler  periodicity.  The 
distance  between  the  spherical  surface  of  observation  and  the 
entrance  to  the  wiggler  is  D.  The  observation  point  is 
located  at  (x,y,z).  Table  I  below  summarizes  some  of  the  FEL 
Interaction  parameters  used. 


TABLE  1. 

Magnetic  wiggler  period  (XJ 

2cm 

Amplitude  of  magnetic  field  (Bw) 

0.05T 

Input  wave  wavelength  (X) 

2Q0Mm 

Input  wave  optical  power  density 

10 7  watts/cm2 

Number  of  magnet  periods 

20 

Electron  energy  (y) 

7 

The  following  restrictions  were  imposed  on  the  calculations: 
a)  no  interaction  between  electrons  were  allowed,  b)  the 
input  wave  field  and  the  wiggler  fields  were  represented  by 
plane  waves,  c)  the  electron  trajectories  were  solely 
determined  by  the  magnetic  wiggler  force,  d)  no  radiation 
reaction  forces  were  included  to  calculate  the  radiation 
fields  produced  by  each  electron  and  e)  no  changes  in 
amplitude  or  phase  of  the  input  signal  were  taken  into 
consideration  (low  gain  conditions). 

All  of  the  results  presented  here  pertain  to  the 
radiation  fields  produced  by  the  electrons  alone.  In  order 
to  calculate  such  things  as  laser  gain  the  radiated  fields 
will  have  to  be  sunned  coherently  to  the  input  field. 


Figure  3  shows  typical  time- records  of  one  of  the 

components  of  the  radiated  field  produced  by  a  beam  of 
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Figure  3«  Time  record  of  the  x-component  of  the  radiated 
electric  field  on  axis  (a)  S  =0Z,  (b)  6  >BZ, 
and  (c)  S z>6j.  2  2  * 

The  dotted  line  represents  the  x-component  of  the 
electric  field  of  the  input  monochromatic  wave. 


electrons  of  length  l=10X  as  observed  on  axis  at  the  point 
(0,0,5m).  The  dotted  lines  represent  the  phase  information 
of  the  Input  wave  (gaussian  TEMoo)  and  the  solid  curves 
represent  the  radiated  fieldl  Figure  3a  shows  the  radiation 
field  observed  for  a  pulse  of  electrons  injected  with  initial 
velocity  equal  to  the  resonance  velocity 
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Initially,  due  to  the  uniform  distribution  of  charge  in  the 
electron  pulse,  there  is  field  cancellation.  As  the  electron 
beam  becomes  bunched  by  the  ponderomotive  force  there  is  a 


finite  growth  of  the  radiated  field.  Very  nearly  at  the  end 
of  the  time  record  the  beam  has  achieved  maximum  bunching  and 
consequently  the  amplitude  of  the  radiated  field  has  reached 
its  maximum  value  also.  Note  that  at  resonance  the  resultant 
radiated  wave  Is  out  of  phase  by  tt/2  with  respect  to  the 
Input  wave.  If  added  to  the  Input  wave  the  radiation  field 
will  not  change  the  amplitude  of  the  Input  wave  (l.e.  no 
signal  gain  occurs).  This  can  be  easily  shown  by  adding  two 
waves  of  the  same  frequency  but  different  amplitude  and 
phase.  If  the  restriction  Is  imposed  that  the  resultant 
amplitude  be  equal  to  the  amplitude  of  one  of  the  Input  waves 
the  following  relation  must  hold: 

tan  ♦  *  ”  2I"  02) 

where  A  Is  the  amplitude  of  both  the  Input  wave  and 
the  resultant  wave,  B  Is  the  amplitude  of  the  second  Input 
wave,  and  is  the  phase  difference  that  must  exist  between 
the  two  input  waves.  For  the  case  in  question  B/A—*>0(low 
gain  condition).  Thus  <fr+-ir/2  as  observed  in  the  time  record 
of  figure  3a. 

The  results  shown  In  figure  3b  represent  the  radiated 
wave  by  the  same  pulse  of  electron  as  in  figure  3a  but  with 
Initial  longitudinal  velocity  @z  >  B*  In  fact,  the  initial 


velocity  was  chosen  to  correspond  with  that  for  maximum 
average  electron  energy  deceleration  (i.e.  at  the  peak  of  the 
one-dimensional  gain  curve).  It  can  be  seen  from  figure  3b 
that  the  radiated  wave  now  has  the  same  phase  as  the  Input 
wave.  Thus,  net  amplification  of  the  Input  wave  Is 
possible.  Similarly,  for  ft  a  net  absorption  of  the  input 
wave  can  occur  (<jwtr).  In  figure  3c  the  radiated  pulse  shown 
was  produced  by  the  same  electron  pulse  but  with  Initial 
0  »8„.  The  phase  of  the  radiated  wave  is  now  oscillating 
in  time.  If  added  to  the  input  wave  the  resultant  wave  will 
have  the  same  amplitude  as  the  input  wave .  The  phase 
oscillation  of  the  radiated  wave  can  be  explained  by  the  fact 
that  although  the  bunching  of  the  electron  beam  occurs  with  a 
spatial  periodicity  close  to  that  of  the  input  wave 

wavelength  the  electrons  radiate  individually  at  shorter 
wavelengths  because  8  Thus,  a  beating  effect  appears. 

To  compare  the  results  obtained  above  with  results 
obtained  In  one  dimensional  FEL  theories,  a  series  of 
computer  runs  were  made  to  determine  the  actual  optical  gain 
curve  on  axis  as  a  function  of  Input  electron  beam  velocity. 
The  radiated  field  was  calculated  on  axis  and  added  to  the 
Input- field  (a  plane  wave)  evaluated  on  axis.  The  on-axis 
gain  was  thus  calculated  and  shown  (see  figure  4)  to  have  the 
same  form  as  that  predicted  by  one  dimensional  theories  and 


GAIN  ON  AXIS 


Figure  4.  Optical  gain  on  axis  as  a  function  of  initial 
longitudinal  velocity. 


that  obtained  experimentally  by  the  Stanford  group.  However, 
one  Important  difference  can  be  noted  Immediately  between  the 
two  results.  The  30  optical  gain  curve  on  axis  Is  shifted  by 
an  amount  53^  toward  larger  Initial  velocity  with  respect  to 
the  ID  gain  curve.  The  30  results  showed  that  the  radiated 
wave  was  In  phase  (♦■0)  with  the  Input  plane  wave  at  the 
point  of  observation.  This  Is  not  In  contradiction  with  the 
results  shown  in  figure  3  because  the  latter  were  obtained 
with  an  Input  gausslan  TEMoo  wave.  The  phase  of  the  TEMoo 
mode  shifts  by  *-90*  far  away  from  Its  waist.  Thus,  at 


resonance  the  radiated  wave  has  a  phase  of  904  with  respect 
to  the  Input  gausslan  wave. 

The  phase  of  the  radiated  wave  is  therefore  very 
important  In  determining  its  coupling  strength  to  the  input 
signal.  Calculations  of  gain  for  input  TEMoo  signals  will  be 
presented  In  a  later  publication. 

ANGULAR  RADIATION  CHARACTERISTICS. 

The  plots  shown  in  figure  5  illustrate  the  angular 
characteristics  of  radiation  in  a  free-el ectron  laser.  Curve 
5a  describes  the  angular  pattern  of  radiation  for  an  electron 
pulse  Injected  with  resonant  velocity  (sz=So)  and  of  very 
short  length  L»x.  The  result  shown  In  the  figure  is  what  is 
expected  from  a  single  radiating  particle.  This  was 
discussed  earlier  (Q^l/y)  and  confirms  the  fact  that  when  no 
coherent  effects  are  present  the  single  particle  radiation 
characteristic  dominate.  The  results  shown  in  figure  5b  were 
obtained  with  an  electron  pulse  of  length  L»50X  and  with 
Initial  resonant  velocity  (S2*6^.  It  shows  how  interference 
effects  reduce  the  solid  angle  Into  which  energy  is  radiated. 
The  half-power  half-angle  shown  is  nearly  35  mrad  compared  to 
60  mrad  for  the  short  pulse.  It  Is  possible  to  think  of  the 
resultant  radiation  as  being  generated  by  a  long  line  of 
phased-array  antenas  radiating  each  short  pulses  of  radiation 
N  periods  long.  From  simple  grating  diffraction  theory  It 
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Figure  5.  Angular  dependence  of  the  radiation  of  a  single 
f l lament,  a)  Very  short  electron  pulse  L«X  at 
resonance  velocity;  b)  Long  electron  pulse  Lm50\ 
at  resonance  velocity;  c)  same  electron  pulse 
injected  with  initial  velocity  0  >6^. 


can  be  shown  that  radiation  will  be  emitted  Into  a  cone  of 
angular  dimension  given  approximately  by  1/No  where  No  is  the 
number  of  radiators  in  the  antena.  For  the  FEL,  however.  No 
represents  the  number  of  electron  bunches  that  can  Interfere 
simultaneously  at  a  fixed  point  In  space.  After  some  though 
It  can  be  shown  that  No*N  for  the  FEL.  That  Is  to  say,  only 
N  optical  periods  can  Interfere  simultaneously  in  the  FEL, 
where  N  is  the  number  of  periods  In  the  wlggler.  Figure  5c 
shows  the  angular  distribution  of  radiation  from  a  beam 
Injected  with  energy  greater  than  resonance  energy.  In  fact 


where  n  Is  the  order  for  constructive  Interference.  Thus,  It 
appears  that  an  electron  beam  prepared  with  the  energy 

necessary  to  give  a  maximum  one-dimenslonally  calculated 

optical  gain  will  In  fact  have  lower  gain  along  the  axis  of 

the  beam.  It  seems  from  figure  5c  that  for  Bz>0*  the 
radiated  energy  best  couples  to  modes  of  higher  order  than 
that  of  the  TEMoo  mode.  In  figure  6c  the  energy  of  the 

electron  pulse  (Ls5Q  X)  was  chosen  to  be  below  the  resonant 

D 

energy  (Sz<3z).  Here  xrad>abunch’  and  as  a  resul-t  there  Is 
no  angle  9  for  which  the  grating  formula  in  equation  (13)  can 
be  satisfied.  Thus,  there  Is  no  ring  structure  in  the 
radiation  cone. 

Figure  7  shows  the  angular  distribution  of  radiation  for 
a  multi  filamentary  beam  injected  into  the  FEL  with  initial 
longitudinal  velocity  sz  =8 z.  One  filament  Is  located  on  axis 
and  10  filaments  were  located  over  a  radius  of  2mm.  All  the 
radiation  waves  were  sumned  coherently.  The  resulting 
angular  radiation  cone  Is  smaller  than  the  ones  shown  in 
figures  5  and  6.  These  additional  interference  effects  are 
not  expected  to  be  seen  in  a  real  experiment  because  of  the 
random  distribution  of  electrons  in  the  beam.  Since  there  is 
no  radial  bunching  force,  the  radiation  fields  emitted  by 
electrons  from  different  transverse  positions  In  the  beam 
should  be  sumned  Incoherently. 


MULTIFILAMENT  ELECTRON  BEAM 
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Angular  dependence  of  the  radiation  from  a  multi- 

filamentary  beam 
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Figure  8.  Angular  dependence  of  the  radiated  wavelength  for 
electron  pulses  of  various  lengths. 


SPECTRAL  RADIATION  CHARACTERISTICS 


The  same  Interference  effects  that  caused  the  PEL  to 
radiate  Into  a  cone  of  smaller  angle  than  that  of  a  single 
electron  are  evident  In  the  radiated  spectrum.  In  figure  8 
the  dependence  of  radiated  wavelength  with  angle  is  shown  for 
beams  of  different  length  injected  with  the  same  initial 
resonant  energy.  It  can  be  clearly  seen  that  as  the  number 
of  radiating  bunches  in  the  beam  increases  the  bandwidth  of 
radiated  energy  decreases.  That  is  to  say,  it  is  possible  to 
put  together  a  string  of  radiated  pulses  and  make  the 
resultant  pulse  more  monochromatic  (i.e.  longer)  than  any 
single  individual  pulse.  Thus,  the  radiation  emitted  by  an 
electron  beam  in  a  FEL  is  monochromatic  by  virtue  of  the 
temporal  interference  effects  that  take  place.  It  is  worth 
noting  that  for  short  electron  pulses  L=A./4  the  radiated 
wavelength  follows  closely  the  curve  for  a  single  electron 
(e.g.  equation  7.) 

CONCLUSIONS 

Results  have  been  presented  that  describe  the  angular 
and  spectral  radiation  characteristics  from  a  constant  period 
FEL  using  the  Llenard-Mlechert  fields.  Interference  effects 
determine  the  angular  and  spatial  characteristics  of  the 
radiated  field.  The  radiated  spectrum  has  been  shown  to  be 


more  monochromatic  and  concentrated  in  a  smaller  radiation 
cone  than  that  for  a  single  electron. 

There  is  sufficient  amplitude  and  phase  Information  in 
the  radiated  field  calculations  to  enable  a  correct 
evaluation  of  the  coupling  of  the  radiated  field  to  the  input 
field.  Results  of  these  calculations  will  be  presented  in  a 
later  paper. 

The  present  computer  code  will  be  shortly  modified  to 
Include  interactions  between  electrons.  In  this  way,  space 
charge  effects  and  radiation  effects  will  be  calculated  in  a 
self-consistent  manner  to  permit  calculations  of 
three-dimensional  radiation  fields  in  the  high  gain  regime. 
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